Constant of motion, Lagrangian and Hamiltonian of the gravitational 
attraction of two bodies with variable mass 



G. Lopez 

Departamento de Fi'sica de la Universidad de Guadalajara 
Apartado Postal 4-137 
44410 Guadalaj Jalisco, Mexico 



PAGS: 03.20. +i 
Octuber, 2005 



ABSTRACT 



The Lagrangian, the Hamiltonian and the constant of motion of the gravita- 
tional attraction of two bodies when one of them has variable mass is con- 
sidered. The relative and center of mass coordinates are not separated, and 
choosing the reference system in the body with much higher mass, it is possible 
to reduce the system of equations to 1-D problem. Then, a constant of motion, 
the Lagrangian, and the Hamiltonian are obtained. The trajectories found in 
the space position- velocity, (x, f), are qualitatively different from those on the 
space position-momentum, (x, p) . 
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1. Introduction 



Mass variable systems has been important since the foundation of the classical 
mechanics and have been relevant in modern physics [1]. Among these type 
of systems one could mentioned: the motion of rockets [2] , the kinetic theory 
of dusty plasmas [3], propagation of electromagnetic waves in dispersive and 
nonlinear media [4], neutrinos mass oscillations [5], black holes formation [6], 
and comets interacting with solar wind [7]. The interest in this last system 
comes from the concern about to determinate correctly the trajectory of the 
comet as its mass is changing. This system belong to the so called two-bodies 
problem. The gravitational two-bodies system is one of the must well known 
systems in classical mechanics [8] and is the system which made a revolution 
in our planetary and cosmological concepts. Normally, one assumes that the 
masses of these two bodies are fixed and unchanged during the dynamical 
interaction [9]. However, this can not be true any more when one consider 
comets as one of the bodies. Comets loose part of their mass as traveling 
around the sun (or other star) due to their interaction with the solar wind 
which blows off particles from their surfaces. In fact, it is possible that the 
comet could disappear as it approaches to the sun [10]. So, one should consider 
the problem of having one body with variable mass during its gravitational 
interaction with other body. 

In this paper, one considers the problem of finding the constant of motion, 
Lagrangian, and Hamiltonian, for the gravitational interaction of two bodies 
when one of them is loosing its mass during the gravitational interaction. The 
mass of one of the bodies is assumed much larger than the mass of the other 
body. Choosing the reference system on big-mass body, the three-dimensional 
two-bodies problem is reduced to a one-dimensional problem. Then, one uses 
the constant of motion approach [11] to find the Lagrangian and the Hamil- 
tonian of the system. A model for the mass variation is given for an explicit 
illustration of form of these quantities. With this model, one shows that the 
trajectories in the space position-velocity (defined by the constant of motion) 
are different than the trajectories on the space position-momentum (defined 
by the Hamiltonian). 



2 



2. Reference system and constant of motion 

Newton's equations of motion for two bodies interacting gravitationally, seen 
from arbitrary inertial reference system, are given by 

d ( dv-i\ Gmim2 . . . 

and 

d f dr2\ Gmim2 , . ,^ . 

where mi and m2 are the masses of the bodies, ri = {xi, yi, Zi) and 

i'2 = {x2,y2, Z2) are the vectors position of the two bodies from our reference 

system, G is the gravitational constant, and 



ri - r2 = r2 - ri 



^J {x2 - xiY + {y2 - yiY + {z2 - ziY 



is the Euchdean distance between the two bodies. It wiU be assumed that mi 
is constant and that m2 varies with respect the time. Taking into consideration 
this mass variation, Eqs. (la) and (lb) are written as 

d^ri Gmima , . 



and 



dh2 Gmim2 , . dv-2. 



where it has been defined m2 as m2 = dm2/dt. Now, let us consider the usual 
relative, r, and center of mass, R, coordinates defined as 

, -o miri+m2r2 

r = r2 — ri , and R — . (4) 

mi + m2 

Let us then differentiate twice these coordinates with respect the time, taking 
into consideration the equations (2) and (3). Thus, the following equations 
are obtained 

(mi + m2)G' 1712 . 

r = r Y2 (5) 

r"* m2 

and 

" -m2 . 2m2m2 . {mi + m2)mim2 - 2mim| 
R = r2 + 7 r + ^ r . (6) 

mi + m2 (mi + m2j (mi + m2) 

One sees that the relative motion does not decouple from the center of mass 
motion. So, these new coordinates are not really useful to deal with mass 
variation systems. In fact, using (4) , one has 

ra = R H ■ r , and = R H ■ r - ■ — r . (7) 

mi + m2 mi + m2 (mi + m2)"' 
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Substituting these expressions in (5) and (6), one can see more clearly this 
coupling, 



nil + m2)G 



+ 



m2{mi + m-i)'^ 



r — 



m2 

1712 



R + 



nil 



mi + ni2 



and 



R 



-111-2 



R + 



/"l/"2 



mi + ni2 (mi + m2)^ 



(nil + 1112)111 iiii2 — niifh^ 
(mi + m2)^ 



(8) 



(9) 



However, one can consider the case for mi ^ m2 (which is the case star-comet), 

— * 

and consider to put our reference system just on the first body (ri = 0). In 
this case, Eq. (3) becomes 



m2 



(Pr Gniini2 



r — m2r , 



where r = r2 = {x, y, z). Using spherical coordinates (r, 9, </?), 
X = r sin ^ cos (p , y — rsmOsinip , z — r cos 9 
Eq. (10) can be written as 



Gniini2 



+ ni2f 



f + 777.2 



(ree 



+ r0 sin ^ ^ ) , 



where f, 9 and (p are unitary directional vectors, 

f = (sin ^cos</7, sin ^sin (/?, cos^) , with f=^^ + 0sin^^. 



and 

(fi — {—sm9,cos9,0) , with (p — sm9 r + cos9 9 . 
Since one has that r = rr, it follows that 

= (r — r^^ + r0 sin^ (/?)f+ (2?^^ + r^ + r(^ sin ^ cos ^)^ 

+ {2rip sm9 + riji sin 9 + 2(p9 cos9)(p , 

and Eq. (12) is discomposed in the following three equations 



m2(r — r9^ + r^sin^ ip) — 



Gmini2 



ni2r , 



m2 (2?^^ -\- r9 -\-r(p sin 9 cos 9) — ni2r9 , 



and 



(10) 



(11) 



(12) 



(13a) 



^ = (cos 6' cos 99, cos 6' sin — sin 6*) , with 9 = —9 r + (p cos9 (p (136) 



(13c) 

(14) 

(15a) 
(156) 



m2 (2r(^ sin 9 + rCp sin 9 + 2ip9 cos 9) = ni2rip sin 9 . (15c) 

Thus, one has obtained coupling among these coordinates due to the term m2- 
Nevertheless, one can restrict oneself to consider the case ni2r ~ 0. For this 
case, it follows that = 0, and the resulting equations are 



4 



(■■ Gmim2 . . , . 

m2[r — rU ) — , [Iba) 

and 

m2{2fe + re)^0. (166) 

Let rrio be the mass of the second body when this one is very far away from 

the first body (when a comet is very far away from the sun, the mass of the 
comet remains constant). Since 1712 7^ on (16b), the expression inside the 
parenthesis must be zero. In addition, one can muhiply this expression by 
nior to get the following constant of motion 

Pe = ruor^e . (17) 

Using this constant of motion in (16a), one obtains the equation 

d'^r _ Gmi ^ P| ma f dr\ 
dt'^ m^r^ m2 \dt J 

This equation represents a dissipative system for m2 > and anti-dissipative 
system for 7712 < 0. Suppose now that 1712 is a function of the distance between 
the first and second body, m2 = m2(r). Therefore, it follows that 



dm2 dm2 dr 
dt dr dt ' 

and Eq. (18) can be written as 



(19) 



^ Gmi ^ mWdrV 
dt'^ m^r^ m2 \dt j ' 

where m'2 = dm2ldr. This equation can be seen as the following autonomous 
dynamical system [12] 

— ^v — = ^ Jj_ - a^y-^ (21) 
dt ' dt m^r^ m2 

A constant of motion for this system is a function K ~ K{r, v) such that the 
following partial differential equation is satisfied [13] 



V 



dK , f-Gmi , P2 m'2 2\dK 



or V m^r"^ 1712 / ov 



This equation can be solved by the characteristic method [14] from which the 
following characteristic curve results 

^/ N 2/ N 2 [ mlir) dr 2P| [ mlir) dr 
C(r, v) = mliry + 2Gm, / -f / ''J , (23) 



7^ fYJj J ^"-j 



and the general solution of (22) is given by 

K{r,v) = F{G{r,v)) , (24) 
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where F is an arbitrary function of the characteristic curve. One can have a 
constant of motion with units of energy by selecting F as F — C/2mo- That 
is, the constant of motion is given by 

K(r v) = !MM^;2 + f ^Ijr) _EL [ "^2(^) dr 

' 2mo ruo J J 



2. Lagrangian and Hamiltonian 

Given the time independent constant of motion (25), the Lagrangian of the 
system (20) can be obtained using the following known expression [11] 

^ / V f Kir, v) dv _ ^, 

L{r,v)^vJ ■ (26) 

Thus, the Lagrangian is given by 

Llr v) = !MM^;2 _ f ^^2(^) dr ^ f ^li^) dr 

' 2mo mo J ml J 

The generalized linear momentum (p = dL/dv) is 

P=^.. (28) 



and the Hamiltonian is 

^_ mop^ , Gmi [ ml{r)dr f ml{r) dr 

2m'.(r)^ J r' mil r^ ' ^ 



2mi(r) 

Note from (25) and (29) that the constant of motion and Hamiltonian can be 
written as 

i^(r,^) = ^-^ + K//(r) (30) 



o 



and 

rUoP^ 



where K// is the effective potential energy defined as 

Gmi f m|(r) dr Pq f m|(r) dr 

This potential energy has an extreme value at the point 



V ^""1 fiT4(j±dr Pi f mi(r) dr 



(326) 



which depends on rUo but it does not depend on the model for m2(r). One can 
see that this extreme value is a minimum for m2(r*) ^ 0, since one has that 

f(fVeff\ _ {GmimoYmoml{r*) 
\ dr^ ]~ P2 ^ 
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On the other hand, because of the expression (28), one could expect different 
behavior of a trajectory in the phase space {r,v) and the phase space {r,p). 
The trajectory r{9) is found using the relation dr/dt = {dr/d9)9, and the Eq. 
(17) in (30) to get 



Je, 



'M=-^r^^M^, (34a) 



where K and Pq are determinate by the initial conditions, K — K{ro, Vo) and 
Pe — f^orlOo- The time of half of cycle of oscillation, T1/2, is obtained directly 
from Eq. (30) as 

fr2 7712 (r) dr 

V^o-Jri - Veff{r) 

where ri and r2 are the two return points deduced as the solution of the 
following equation 

Veff{n)^K, i = l,2. (34c) 



3. Model of Variable Mass 



As a possible application of (25) and (29), consider that a comet looses material 
as a result of the interaction with star wind in the following way (for one cycle 
of oscillation) 



m2[r) 



niooV^ — e incoming {v < 0) 

^ mjc"^''^"''^ + m/(l — e""'") outgoing (v > 0) 



(35) 



where rrioo or m/ (where mj = 2mj — nioo by symmetry) is the mass of the 
comet very far away from the star (in each case), nii is the mass of the comet 
at the closets approach to the star (distance ri), = mooVl — e~°''i, and a 
is a factor that can be adjusted from experimental data. Thus, the effective 
potential (32a) has the following form for the incoming case (mo = moo) 



(in) 



r) — 



Gmirric 



-(1-e- 



+ 



P2 



^runnr^ 



:i-e- 



+ 



GMiniooOt + 



a^Pl 



2m„ 



Ei{—ar) + 



aPie 



2mo„r 



(36a) 



where Ei{x) is the exponential-integral function [15]. 
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For the outgoing case, one has nio — ruf and 



{out) 



eff 



(r) = 



Gminif ^ Pq 



ruf 



-2ar 



m 



f 



r 

-2ar 



2aEi{-2ar) 



2r2 



+ 



ae 



-2ar 



+ 2Ei{-2ar) 



+2Gmi(mie""^ -m/) 



ml 



r 

-ar 



— aEi{—ar) 



2r2 



+ 



ae 



a 



2r 



+ —Ei[—ar) 



(366) 



where is defined now as Pe — mfr'^9. The extreme point of the effective 
potential (32b) for the incoming and outgoing cases is given by 



(37) 



Given the definition (35), the constant of motion, Lagrangian, generahzed 
finear momentum, and Hamiltonian are given by 



2m„ 



and 



m|(r) 



2mi(r) ^ ^ff^^ ' 



(38) 

(39) 
(40) 

(41) 



where i — in for the incoming case, and i — out for the outgoing case. As 
an example of illustration of this model, let us use the following parameters 
to estimate the dependence of several physical quantities with respect the 
parameter a, 

G = 6.67 X IQ-^^m^/Kg sec ; rUoo = 10^ Kg ; m/ = 0.1 rUoo ] 
Pe = lO^^Kg rri^/sec ; and X = -8 x lO'^^ Joules . (42) 

Fig. 1 shows the curves of Veff{r) for several values of a (incoming case). As 
one can see from this figure, the location of the minimum does not change, 
but the minimum value of V^ff tends to disappear as a goes to zero. Also for 
the incoming case. Fig. 2 shows how the minimum distance of approximation 
of the two bodies, ri, and maximum distance, r2, behave as a function of the 



8 



parameter a. As one can guess, the following limit is satisfied lim^^o^"! = 
lima^o^2 = which will become a inflexion point for Vg//. Fig. 3 shows 
the velocity [v] and normalized linear momentum {p/mo) as a function of r 
for several values of a and for the incoming case. All the trajectories start at 
r2 = 200 and finish at ri{a). One can see the difference of the trajectories in 
(a) with respect to (b) due to position dependence of the momentum, relation 
(40). 

Conclusions 

The Lagrangian, Hamiltonian and a constant of motion of the gravitational 
attraction of two bodies when one of them has variable mass were given. One 
found that the relative and center of mass coordinates are coupled due to 
this mass variation. However, chosen the reference system in the much more 
massive body, it was possible to reduce the system to 1-D problem. Then, the 
constant of motion, Lagrangian and Hamiltonian were obtained. One main 
feature of these quantities was the appearance of an effective potential, which 
is reduced (when m2 = 0) to the usual gravitational effective potential of two 
bodies with fixed masses. Other feature was the distance dependence of the 
generalized linear momentum. A model for comet-mass-variation was given 
which depends on the parameter a. A study was made of the dependence with 
respect to a of V^//, minimum and maximum distance between the two bodies, 
and the trajectories in the spaces {r,v) and (r,p). Of course, the problem of 
the interaction comet-star with the variation of mass deserves more complete 
analysis. The intention here with this example was to show explicitly the form 
of the constant of motion, Lagrangian, and Hamiltonian and to point out the 
different trajectories behavior in the spaces (r, and (r,p) arising from the 
constant of motion and Hamiltonian. 
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Figure Captions 

Fig. 1 v}jl\r) with the values of the parameters given on (42), for a = 1 (1); 
a = 0.01 (2); and a = 0.005 (3). 

Fig. 2 Maximum (r2) and minimum (ri) distances between the two bodies as 
a function of the parameter a. 

Fig. 3 (a): Trajectories in the plane (r, f); (b): Trajectories in the plane (r,p). 
a = 1 (1), a = 0.01 (2), and a = 0.005 (3). 
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